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THE  EFFECT  OF  THICKNESS  OH  AIRFOILS  WITH  CONST.iNT  VERTIC.X 
ACCELERATION  AT  SUPERSONIC  SPEEPS 


ABSTRACT 

The  effects  of  thickness  on  the  lift  and  pitching* moment  of  two 
dimensional  airfoils  at  supersonic  speeds  with  constant  vertical 
acceleration  are  investigated.  The  airfoils  considered  have  arbitrary 
syrmnetrical  cross  sections,  and  the  flow  is  supersonic  everywhere. 


The  analysis  is  based  on  a second  on_er  theory  similar  to  the 
second  order  theory  introduced  by  Buseraann  and  extended  by  Tan  Dyke. 
The  lifting  pressure  due  to  a constant  vertical  acceleration  is  found, 
and  this  is  used  to  calculate  the  lift  coefficient,  CT  , and  the  momeflt 


coefficient,  C_  , 
* m* 
a 


V 

due  to  a constant  vertical  acceleration. 


* The  airfoil7 s second  order  contribution  to  the  damping  of  longitudi- 
nal oscillations  in  aircraft  is  considered  and  its  relation  to  the 
damping  of  oscillating  airfoilfr  is  investigated.  • 
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V 


INTRODUCTION 


The  development  of  the  linearized  theory  of  supersonic  flow  has 
permitted  a first  order  evaluation  of  a number  of  stability  derivatives. 
Second  order  theories  similar  to  the  one  introduced  by  Busemann  (ref.  1) 
and  extended  by  Van  Dyke  (refs.  2 and  3)  offer  possibilities  of  obtain- 
ing second  order  evaluations  of  certain  stability  derivatives,  such  as 
the  lift,  CT  „ and  the  moment.  C , due  to  a constant  vertical  accelera- 
Lc  ma 

tion.  Examples  of  the  use  of  second  order  theories  to  evaluate  stability 
derivatives  are  found  in  refs.  (U)  and  ($). 

Recently  two  pap>er3  have  appeared  (refs.  (6)  and  (?))  dealing  with 
two  dimensional  oscillating  airfoils.  The  flow  over  oscillating  airfoils 
is  of  great  importance  in  flutter  calculations.  For  stability  studies, 
however,  it  is  more  convenient  to  knew  the  flow  cr  jr  airfoils  with  a 
constant  vertical  acceleration  and  a constant  rate  of  pitch. 


In  thi3  paper  a second  order  theory  is  developed  for  two  dimensional 
airfoils  with  constant  vertical  acceleration  at  supersonic  speeds. 

From  this  theory  the  lifting  pressure  due  to  constant  vertical  accelera- 
tion is  evaluated.  The  expression  for  the  lifting  pressure  enables  the 
stability  derivatives  and  ^ to  be  calculated.  The  airfoil3  con- 

a a 


3idered  herein  have  arbitrary  symnetrical  cross  sections;  however,  the 
analysis  can  easily  be  extended  to  include  airfoils  with  unsynmetrical 
cross  sections. 


The  flow  around  an  accelerating  airfoil  is  of  an  unsteady  nature 
and  this  necessitates  the  U3e  of  time  dependent  partial  differential 
equations.  The  differential  equation,  studied  here  Is  not  the  time 
dependent  equation  coranonly  used,  since  the  coordinate  axes  employed 
in  the  analysis  are  attached  to  the  airfoil. 

analyse: 


Introduction*  Recent  work  by  Milton  D.  Van  Ityke  (ref3.  2 and  3) 
indicates  that  seOond  order  s' ^utions  of  the  partial  differential 
equation  of  steady  supersonic  flow  can  be  obtained  by  iterative  methods. 
In  the  present  paper  it  will  be  assumed  that  iterative  methods  can  also 
be  used  to  obtain  second  order  solutions  of  the  differential  equation 
for  unsteady  supersonic  flow.  This  assumption  wa3  P"«de  in  ref.  (6) 
in  treating  oscillating  airfoils. 

In  addition  it  will  be  assumed  that  the  characteristics  are  the 
same  for  the  first  and  second  order  solutions*  This  assumption  does 
not  appear  unreasonable  3ince  for  steady  plane  flow  the  second  order 
solution  (ref.  (2))  found  by  an  iterative  method  yields  the  correct 
second  order  pressure  of  the  Busemann  second  order  theory. 
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The  Partial  Differential  Equation:  The  partial  differential  equation 

for  the  second  order  potential  flow  around  an  airfoil  is  (from  ref.  (U)) 

•>***•«•  4*  ,•«  • \ f t t * N2f(r-i)($x..  * 


?<£  5: 

^xf  ^x'x* 


2 4>  5 

z - 


X 2 


V r X 


X t 


\T  ^z’  *z't  1 


(1) 


The  perturbation  potential  function,  'J-  « has  been  normalized  through 
division  by  the  free  stream  velocity  Also  from  ref  (U)  the  pressure 
relation  is 


Cp  * ' 2 (^x  /V)  *1  * x ^ * 21(2 ?x-  $t’/V  * M^t  (2) 


These  two  equations  are  associated  with  axes  fixed  in  space*  For 
the  problem  considered  here  it  is  convenient  to  express  the  differ- 
ential equation  in  terns  of  « new  set  of  axes  fixed  to  the  airfoil. 

The  relations  between  the  two  sets  of  axes  are  (see  fig.  1) 


x * x 
Z - 2 

t * t 


aV  (t  )2 


(3) 


For  the  axes  attached  to  the  airfoil  eqs.  (l)  and  (2)  become 

-(2M2$xt/V)  (^ft/c2)  - 2M2|it#„  ♦ 


xt' 


tr 


xz 


C - 


(at4zt/V)  ♦ j «§Z/(2V)J  ♦ Uy-i)/2  ($x  /V)  ($^4  $zs)  ♦ 

$ $ *£  & yk  ♦/fi  « ,/v\l  (U) 

x^xx  *z  *xz  , x *"xt  j s~Z  * zt'  j * v ' 

o 

-2jsx  /?)♦  *P2f/  * ( 2X2 <JX  §t/v ) * (5) 

(«2#t2  /7s) 


Bqs,  (U)  and  (5)  will  be  used  m the  analysis  of  accelerating  airfoils. 

Solution  By  Iteration  On  the  assumption  that  eq»  (U)  can  be 
solved  by  an  iterative  procedure  the  initial  step  will  be  finding  the 
first  order  solution.  The  first  order  partial  differential  equation, 
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obtained  by  neglecting  the  second  order  terms  in  eq.  (d),  is 

- P2  <3  ♦ 0 - (2M2  0 ./V)  (0ftA2)  * 0 (6) 

XX  ZZ  XX  tt 


The  solution  to  eq.  (6)  is  taken  as  the  first  approximation  to  the 
solution  of  eq.  (U)  We  shall  assume  that  the  second  approximation 
can  be  found  by  substituting  the  first  order  solution  into  the  right 
side  of  eq.  (h)  and  solving  the  resulting  non- homogeneous  equation, 
which  from  eqs.  (U)  and  (6),  is 


- 2m2{«  dxi  .(it  0jt/V). 

[i  V(JV,j  * [<r  ' a/x*‘ 

ax  **xx  * ' "‘z 


(7) 


The  solution  of  eq  (7)  will  be  referred  to  as  the  second  order  solution. 


Examination  of  Form  of  Potential  Function:  It  is  helpful  to  in- 

vestigate  the  type  of  solution  obtained  from  eq.  (7)  The  first  order 
solution  will  be  cf  tne  form  0 * a g,(x  i)  ♦ a g (x.z.t)  ♦«g,(x>*)  ? 
where  t.  is  a thickness  parameter  and  &iere  o.  a and  c are  small  compared 
tc  unity.  It  follows  from  the  above  expression  and  eq  (7)  that  the 
second  order  solution  wiU  be  of  the  form 

V-  a^h^vx  z^  * aa  h^'x  z t)  ♦ a'  h^(x  z ,t)  * afh^x  z)  + 
d^h^fx  2 t)  hfo(x  ?) 

This  paper  Is  concerned  only  with  the  lifting  pressure  and  its 

integrated  effects.  The  terms  h-j(x.z),  a ah^Cx  z.t),  and  o^h^(x,z,  t) 

contribute  nothing  to  the  lifting  pressure  This  can  be  seen  from  the 
following  argument  The  thickness  parameter  t is  not  present, 
therefore,  the  airfoils  can  be  considered  to  have  zero  thickness  when 
these  terms  alone  are  treated.  If  the  potential  cf  the  flow  on  the 
upper  surface  of  a flat  plate  is  expressed  as 

$(a.a)  -1-  og1  r dg.,  - a^h,  * aah^,  ♦ o^h^ 

the  potential  or.  the  lower  surface  is  giver,  by 

* / » x • 2 • * 2 
9 ( a -a)  - o g,  a g?  * a h^  ♦ aah?  ♦ a h^ 

The  potential  difference  is 
« 2agT  ♦ 2dg^ 
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Since  for  the  flat  plate  the  pressure  difference  between  tne  upper  and 
lower  surfaces  can  be  found  directly  from  the  potential  difference,  the 


2 • - 2 

terms  a a o and  a do  not  affect  the  lifting  pressure. 

The  terms  on  the  T*ight  side  of  c-q.  (7)  which  are  multiplied  by 

2 • • 2 /\  2 
a t an,  and  a can  be  reflected;  fcr  «q.  (.7)  is  linear,  and  a h^.aah^, 


m 2 

and  a h.  do  not  contribute  to  the  lifting  pressure.  The  terms 

* 2 » 

aeh^  and  £ h^  are  independent,  cf  a ( the  acceleration  parameter), 

and  hence  they  also  contribute  nothing  to  the  pressure  due  to 

„ • 

acceleration 


The  remaining  term,  art  ' x i t)  will  be  found  by  use  of  eq,  (7)y 

> 2 2. 

neglecting  all  expressions  involving  a a-'  . a and  aa 

A further  consideration  of  the  form  of  the  solution  indicates 
that  the  second  order  lifting  pressure  is  linear  m the  thickness 
parameter.  Tnis  can  be  established  by  considering  a first  order 
solution  of  the  fcrm 

Q * a g^(x  z)  * ag^(x  z t)  • < x 2)  *■  r ^ g^(x  1) 

Then  the  second  order  lifting  pressure  will  be  cf  the  form 

C ’ I (xt)  ♦ a c.  T,  k t ) * a r * . ( x : ) * ... 

This  equation  is  linear  1 r.  r anl  thus  the  lifting  pressure  for 

various  known  thickness  d:  stnou  tiers  lan  be  added  tc  obtain  the 
lifting  pressure  for  other  thickness  distributions. 

This  linearity  has  no  spe^ia.  value  fc’  the  *vr.  dimensional 
airfoil,  sirte  we  shall  de’ermir-e  the  solution  for  an  arb.trary  thickness 
distribution.  It  can  be  shown  similarly  that  linearity  of  lifting 
pressure  with  thickness  du?.  r ifci*  ;c:  hold;  frr  three  dimensional  air 
foils  also,  and  for  this  case  it  shcu.d  prove  quite  usefxu. 

Boundary  Cond’firr.s  Phvsica.  considerations  require  that  the 
flow  be  tar  gen  tc.  the  surface  cf  ’he  airfoil  ar.d  tliat  all  velocity 
perturbations  vanish  upstream.  *.f  tr.e  airfoil.  These  boundary  con- 
ditions may  be  expressed  mathemat; cally  as 

0 (x  y z t ) - 0 upstream  cf  the  airfoil. 

'••  ( x y 2 t ) • 0 

and  v "7^  * 0 

* The  terms  a >'  < u h^  and  a4"  h are  asso  lat.ed  with  steady  super- 

sonic flew  and  c.  ;r.  fc?  found  by  she  Bj remain  second  order  theory. 


where  s(x^  a)  "0  13  the  equation  of  the  jut  face  of  the  atrfc.  1 


The  equation  of  the  surface  of  the  airfoil  nay  also  be  expressed 
as  z * < f(x) 

Thus 

v ■ • uf' 

s 

Since  the  velocity  vector  v.  can  be  written  as 


v x 1 (l  * 0 ■*  x.  ) ♦ r ♦ v ♦ at) 

xx  z z 

it  follows  that  the  boundary  condition  cn  the  airfoil  surf  ce  is 


(1  ♦ 0 ♦ ) f.  f * 0 - ♦ at  * 0 

xx  z z 

The  coordinate  axes  wi.i  be  chcser  so  th.it  the  airfoil  lies 
approximately  in  the  \ - 0 plane..  The  boundary  conditions  for  the 
first  order  solution  are  0(x  y 2)  C upstream  of  the  arfoil  and 


0 

7 


* a 


at  - f 


(8) 


Since  only  the  aC  term  is  oeing  calcu.ated  a wi  1 1 be  set  equal  to 
zero.  Similarly  th«*  boundary'  '-onditicns  fcr  the  second  order  solution 
are  \p  - 0 upstream  of  the  airfoil  and 


2 - 0 


< f 


Q> 


1 1 


(9) 


For  the  airfoils  considered  m this  paper  the  first  order  velocity 
components  are  discontinuous  across  the  Macn  shee'  from  the  leading 
edge.  The  effect  of  these  discontinuities  on  the  se*.or.d  order  solution 
will  be  determined  by  assuming  that  a small  transition  surface  is 
attached  to  the  leading  edge.,  Tnio  small  surface  is  sc.  shaped  that  the 
discontinuities  in  the  velocity  components  are  reroevea.  The  effect  of 
the  discontinuities  will  be  evaluated  by  a limiting  process  in  which 
the  width  of  the  attached  surface  approaches  zeio.  A surface  which 
will  remove  tne  first  order  velocity  discontinuities  in  the  flow  over 
the  upper  surface  is 

1 ' - (a  x^t /2)  * ( x*  a. /«?) 

where  the  surface  extends  from  x - 0 to  x • t , and  a/~a>a/*c 

Evaluating  the  se.or.d  order  potential  function  associated  with  the 
flow  over  the  snlaii  transition  surface  introduced  above  will  lead  to 
the  determination  cf  the  second  order  loading  edge  di see ntinui ties 
The  first  order  solution  fcr  the  f.cw  ever  the  transition  surface  is 

r ■>  If',  n 0 1 

0 * j (x  pz)  /(?pv) | ‘ ' a r(x  U /<  3p")  * oVt  • Va  ; 


(ID) 


The  second  order  solution  can  be  found  by  substituting  eq.  (ID) 
into  eq.  (7)  and  solving  the  reulting  non-homogeneous  equation; 
however,  in  thi3  case  it  is  easier  to  use  a different  approach.  In 
appendix  A it  is  shown  that  the  solution  for  planar  problems  to  the 
three  dimensional  non-homogeneous  equation 

-P2^  + Yyy+VZz~  (2M2yxt/V)-  (*tt/c?)  - F(x,y,*,t)  (ll) 
can  be  written  as 

y(x,y,*,t)  • <J  S ^ !_r(ti)  ♦ F(t2)"j  dv  ~ (12) 

o ~ 

1 ff  V- 

* 0 (. 

where  vQ  is  the  volume  inclosed  by  the  forward  Kach  cone  from  the  point 
(x,  y,  z) , and  the  surface  Sq  is  the  area  of  the  transition  surface  in 
the  forward  Mach  cone  from  the  point  (x  y,  2) . 


When  eq.  (12)  is  applied  to  an  arbitrary  point  located  on  the 
Mach  surface  from  the  trailing  edge  of  the  transition  surface,  the 
result  can  be  expressed  as 


V^(x,*,t)  - w 


,SJ 


p:  ♦ ( 


{J(x-£)2- p2(z-{)2)/p 


PS 


dr)- 


l_ 

2n 


r( 


d4 


-(Ax-l)2  -p2U-£)2)/p 
(\j(x  -£)2  -P2z2  )/p 

-t  /Rra5/N'00'ti,M 

7$  r 


d,| 


Sine.  F,  [dy(5.  0,  0,  tj)  /a{],  and  [aft  0,  0,  t2)/aC] 
independent  of  Y)  , the  preceding  equation  can  be  reduced  to 

vu.z.t)  - \ r1  ic  rPJ'f  [«$, 

o 13 

j j {[n<l  • o o-V/  «{]♦  [at<£o,<M^»([jd$ 


j 


y 


12 


as 


The  integrand  of  the  first  integral  in  eq.  (13)  can  be  written 
*(£.,(>  V ♦ t2)  - J0(5.$.  tW£,-p^)Ji(4,^,t)  ♦ 

(4-  K)2  J2<C.C.  o ♦ (4-  fK>3  j3<4,£,  t) 

are 

where  JQ,  J^,  «I>  ?and  .continuous  functions  and  are  not  zero  at£"  p£* 

Since  the  4 in  the  integrands  of  eq.  (13)  varies  over  a very  small 
range,  the  functions  Jq,  J^,  and  can  be  considered  as  constants 

with  respect  to  integration  in  the  4 direction.  Therefore  ^ can  be 
replaced  by  p £ in  these  integrands,  thus  permitting  integration  to  be 
performed  with  respect  to  4 • Eq.  (13)  becomes 

y(x,*,t) / <ypC>£.-t)  J jx(p^, t)  d$  ♦ 

0 o 

■fof  Vo  c . »>  * ilj  f j3ik  . c - 1>  dr  <^> 

.2  f d M2  / X / ^ O’  Vt  a 1 

~r  *!J 

The  value  of  the  second  order  potential  function  on  the  downstream 
side  of  the  leading  edge  Mach  sheet  above  the  airfoil  can  be  found  by 
taking  the  limit  of  eq.  ( UU)  as  / approaches  zero  and  as  d and  a 
approach  infinity  so  that 

lin  of  • & 

/-o 

a — ► °o 


and 


lim  % / ■ £ 


? 

a 


-o 
. <ao 


Thus  the  discontinuity  in  the  second  order  potential  function  is  given 
by 

^x-pz  f+o  l£~  J J0(p£>  t)  d£*  ^ / J1(P<,?,t)d^  (i5) 


a-*»oo 


c ■♦•oo 
a -►00 


The  integrals  containing  and  are  zero  in  the  limit  since  Jg  and 
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2 • 2*2  * 

are  linear  functions  of  a , aa,  a , a , and  aa . 

From  eqs.  (ID)  and  (7)  the  functions  J and  J,  are  (neglecting 
2*2*2  ox 

the  bc , a a,  a and  a*  terns) 

JQ  - l*d  M2  a [ -t  ♦ M(x-p  £ ) / (p2c)] 

[vt-  hVp2]  /(vp2) 

Substituting  the  preceding  expression  into  eq.  (15)  and  performing 
the  indicated  operations  yield 

u/|  - «*(M2/p2)  f-t  ♦ K x/(2P2c)1  pz  ♦ 

r X - p2  L J 

a £ [(y^D/^m^/P1*)  [t-dx/(p2c)]  pc  (16) 

Expression  (16)  is  the  value  of  the  discontinuity  in  the  second  order 
potential  function  across  the  leading  edge  Mach  sheet  above  the  airfoil. 

Solution  of  the  Partial  Differential  Equation;  The  part  of  the 
second  order  potential  function  which  contributes  to  the  lifting 
pressure  due  to  a constant  acceleration  will  now  be  determined.  Since 
the  method  of  solution  is  essentially  the  same  for  both  the  upper  and 
lower  surfaces,  only  the  flow  over  the  upper  surface  will  be  considered 
in  detail. 


The  first  order  solution  is 

0 - ♦ y * * '€  (17) 

It  follows  from  eqs.  (7)  and  (17)  that  the  second  order  potential 
function  must  satisfy  the  non-homogeneous  equation 


- f*xx"**z 


2M 


xt 


1 M2  0€ 

”2  ^tt  " T 

c 


r 

.i  ♦ 


yM^i-p2) 


qU 


f' 


[(-fM2  - M2  ♦ 2)/p2  ] Vtf"  ♦ (H?/pli)(r+l)(x*p3  i)  f"j  (18) 
2 • 2 *2 

where  the  a , a a,  a e , c , and  o terms  have  been  neglected. 


The  potential  function  will  be  divided  into  two  parts  and  0^ 
such  that  0^  satisfies  the  non-homogeneous  equation  and  6^  satisfies 
the  homogeneous  equation. 


11* 


By  Inspection  a solution  of  the  non-honogeneous  equation  is  found 

to  be 

e - ■ £ ISL  (M2  y - M2  ♦ 2)  ztf * - ££i£  [Vo1*  - 2)4H**  - 3*1  *f  - 
1 2pJ  Uv  y L J 

[pk  - Y I f ♦ (r+D  xzf‘  (19) 

UV  pU  L J 2V  p5 

The  boundary  conditions  are  given  by  eqs.  (9)  and  (16).  From  eq. 

( 19)  and  the  boundary  conditions  it  follows  that 
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By  inspection 
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f(u)  du 


} 


(20) 


Thus,  from  eqs.  (19)  and  (20),  the  part  of  the  second  order 
potential  function  which  contributes  to  the  lifting  pressure  is 

Y - — j fn2(n2y  - M2  ♦ 2)  p z Vtf  - y ^ P2  *2  f -(M^/p2V^i)3tpf'z  + 
2Vpu  t 'V  / 

(h\-  - ♦ 2)  Vtf  -(M2/p2)fclAY  - ♦ M2  ♦ l)pz  ♦ ('m2y-M2+2)  xl  f - 

t - fir.  •%  « 


(V/P2)^  Y " 3H1*  ♦ -5M2  - l)^  Pz  f(u)  du  | 


(21) 


as 


Lifting  Pressure: 


The  lifting  pressure  distribution  can  be  expressed 
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From  eqs.  (2)  and  (21)  it  follows  that  for 
lifting  pressure  distribution  at  t * 0 is 
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o Ux 
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V fT 
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♦ 2) 


symmetric .1  airfoils  the 


x f»(x)  ♦ (2M2y  -3M2+5)f(x)] 

(22) 


Fig.  2 presents  the  lifting  pressure  distribution  at  t*0  on  a 
ten  percent  thick  wedge  for  various  Mach  numbers,  and  fig.  3 presents 
the  lifting  pressure  distribution  at  t ■ 0 on  a five  percent  thick 
airfoil  with  a parabolic  cross  section  for  various  Mach  numbers. 


The  Force  and  Moment:  The  effect  of  thickness  on  the  lift, 

and  the  moment  C , due  to  a constant  vertical  acceleration  can  be 

“d 

found  by  use  of  eq.  (22).  The  stability  derivatives  and  Cm> 

a a 

can  be  expressed  as 


From  the  preceding  relations  and  eq.  (22)  the  C.  and  the 


C of  an  airfoil  are  found  to  be 
m • 


h ' F 
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Fig.  U presents 

with  Mach  number  for 
variation  of  the  C. 

Ld 

Mach  number  for  xero 
the  variation  of  the 


the  variation  of  the  CT  of  a wedge  airfoil 

L& 

various  thicknesses^  and  fig.  5 presents  the 
of  an  airfoil  of  parabolic  cross  section  with 

and  7.5  percent  thicknesses.  Fig.  6 presents 
C of  a wedge  airfoil  with  Mach  number  for 
nd 


various  thicknesses,  and  fig.  7 presents  the  variation  of  the  C 

ma 


of  an  airfoil  of  parabolic  cross  section  with  Mach  number  for 
various  thicknesses. 


The  sum 


of  the  stability  derivatives  C 

TO  a 

d 


and  the  damping  in 


pitch,  Cra  , largely  determines  the  damping  of  the  longitudinal 

oscillations  in  aircraft.  The  second  order  damping  in  pitch  far 
airfoils  with  sjmnetrical  thickness  distributi jns  is  given  by 
eq.  (23)  of  ref.  (^).  The  sum  of  this  equation  and  eq.  (2U)  can 
be  written  as 


(d/cr)(*W-WJ*u)](i-  i)  [WW-W 

_ C " 

18M2  - 8jJ  dx  ♦ 2(mS^-UM2+U)  J^1*  xf(x)  I 


11M^  ♦ 

(25) 


Fig.  8 presents  the  variation  of  the  cn.  * ^ of  a wedge  airfoil 

a q 

with  Mach  number  for  various  thicknesses  for  the  case  where  the  axis 
of  pitch  is  located  at  the  mid-chord  poim,.  Fig.  9 presents  the 
variation  of  the  C ♦ C of  an  airfoil  with  a parabolic  cross  section 

W*  TO 

d q 

with  Mach  number  for  various  thicknesses  for  the  case  where  the  axis  of 
pitch  is  located  at  the  midchord  point.  These  two  figures  indicate  that 
the  effect  of  thickness  has  a destabilizing  effect  for  a wedge  airfoil 
and  has  a stabilizing  effect  for  an  airfoil  with  parabolic  cross  section. 


Fig.  10  presents  the  regions  of  possible  instability  for  a ten 
percent  wedge  and  a five  percent  parabolic  arc  airfoil.  The  curves 
are  lines  of  zero  damping  found  by  placing  C ♦ Cm  equal  to  zero. 

a q 

This  figure  indicates  that  the  effect  of  thickness  increases  the  region 
of  instability  for  the  wedge  airfoil  and  shifts  the  region  to  a lower 
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value  of  d/cf  for  the  parabolic  airfoil. 


The  relation  between  C ♦ C 

nq 


and  the  damping  of  a slowly  oscil- 


lating airfoil  is  investigated  in  appendix  B.  It  is  shown  that  to 
the  second  order  in  amplitude  and  thickness,  and  to  the  first  order  in 
frequency,  the  damping  of  a slowly  oscillating  airfoil  is  directly  pro- 
portional to  the  sura  C ♦ C 

ra. 


m_ 


CONCLUDING  REMARKS 

The  airfoils  considered  in  this  paper  have  syraraetrical  thickness 
distributions.  But  since  the  flow  over  the  upper  and  lower  surfaces 
of  the  airfoils  treated  are  independent  of  each  other,  the  aerodynamic 
properties  due  to  constant  vertical  accelerations  of  airfoils  with 
unsymnetrical  thickness  distributions  can  easily  be  determined  from 
the  results  obtained  here. 

The  limitations  of  the  Buseoann  second  order  theory  have  been 
investigated  (see  ref.  8).  Since  the  theory  contained  in  the  present 
paper  is  closely  associated  with  the  Busemann  second  order  theory 
it  seems  likely  that  the  results  presented  herein  have  similar 
limitations. 


Qo*k/v\,  C. 
John  c.  martin 
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APPENDIX  A 


In  ref.  (9)  several  expressions  are  given  for  the  scalar  potential 
function  of  the  three  dimensional  time  dependent  linearized  partial 
differential  equation  of  supersonic  flow.  In  this  appendix  the  results 
of  ref.  (9)  are  extended  to  include  eq.  (11).  The  notation  of  ref.  (9) 
will  be  used  in  this  appendix. 

Eq.  (17)  of  ref.  (9)  can  be  written  as 


where  the  closed  surface  of  the  integral  of  the  left  side  incloses  the 
volume  of  the  integral  on  the  right  side. 


The  potential  function,  0,  is  required  to  satisfy  eq.  (ll).  In  this 
case  eq.  (A-l)  reduces  to 
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Eq.  (A-2)  will  be  applied  to  a volume  (denoted  by  v ) inclosed  in 
the  forward  Mach  cone  from  the  point  (x,y, z)  (see  fig.  (ll)).  This 
volume  is  bounded  by  the  forward  Mach  cone,  an  arbitrary  surface, 

inclosed  in  the  forward  Mach  cone,  and  a surface  given  by  ^ ■ x - 6, 
where  6 is  small. 


The  finite  part  of  the  surface  integral  is  zero  since  on  the 
Mach  cone 


tl  " *2 

The  surface  integral  over  the  area 

£ ■ x -6 

reduces  to  a time  independent  problem  in  the  limit  as  6 approaches  zero. 
Thus  this  surface  integral  in  the  limit  becomes 
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The  details  of  the  above  integration  aa y be  found  in  ref.  (10)  or 
ref.  (11). 


It  follows  from  the  preceding  paragraph  and  eq.  (1-2)  that  tbs 
potential  function  at  the  point  (x,y, *)  can  be  expressed  as 

0(x,y,«,t)  - - ji-  [(Fj  ♦ F2)/r]  d.  * 7J- VI- 


U-3) 


For  planar  probleas  in  which  the  pqint  (x,y,^) 
plane  and  the  disturbing  surface  lies  in  the  t 
can  be  reduced  to 

j*(x,y,z,t)  - - X—  JI^F1  + F2)/R]dT  " 75“ 
R ^1  R ^1  R-  ^2  R 

?;  ^ 7^  71  *i'7l 


lies  above  the  z ■ 0 
* 0 plane  eq.  (A-3) 


^-)]df,  d>, 


where  the  surface  is  the  z ■ 0 plane. 

The  potential  function  is  as  7et  undefined.  The  potential 
function  will  be  defined  so  that 

0'(x,y,-lo,t)  - ft(x,y,&o, t) 

In  this  case  eq.  (A~U)  reduces  to 

0(x,r,.,t)  - - i jJJ  [<  V V*]  d.  - < 

ve  S,  s 


d7j 

(A-5) 
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APPENDIX  B 

The  Relation  Between  the  C Plus  the  C and  the  Damping  of 

“4  mq 

Oscillating  Airfoils 

The  damping  of  oscillating  airfoils  calculated  on  the  basis  of  the 

linearised  theory  can  easily  be  shown  to  be  directly  proportional  to 

the  sum  C ♦ C calculated  by  the  linearized  theory.  The  question 
m.  m 

o q 

arises  — does  this  hold  true  for  second  order  calculations?  This 
question  will  be  answered  by  investigating  the  partial  differential 
equations  apd  the  boundary  conditions  for  the  three  types  of  motions 
involved. 

Eq.  (1)  can  be  used  in  the  analysis  of  oscillating  airfoils;  how- 
ever, for  comparison  purposes  it  is  more  convenient  to  express  the 
equation  in  terns  of  axes  fixed  to  the  airfoil.  The  relations  between 
the  two  sets  of  axes  are  (see  fig.  12) 

x ■ x'  cos  e ■*or*’  - z1  sin  e 

z ■ x'  sin  a e'^'  ♦ z1  cos  o (B-l) 

m m 

t - t' 

For  the  axes  attached  tr  the  airfoil  eq.  (1)  becomes  (to  the 
second  order) 

♦$„  - -(l/=2)  *tt  - o.  e [»*[#„*(** 

(iun.  fjo/V)  x^/V)  ♦ (i »§Z/V 

(**/<2)(i  ?xt  * x $ zt)  4 (“2/c2J(*  $ X " *^z)]*  (&"2) 

K?  L(irl)fe  * St  /v)(<P„  *i„)  * ♦ 2$,  *„  * 

(JM*t/r) 

For  slowly  oscillating  airfoils  the  time  variation  can  be  ex- 
panded in  the  form 

e i*t  - 1 ♦ i*t  - »2t2/2  ♦ . . . 
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The  result  of  substituting  this  relation  into  eq.  (B  2)  and  retaining 
only  the  terms  up  to  first  power  of  <o  yields- 

-e2  $ » ♦ z *,/'*)■  **tt/v)* 

am  iorfM2  ft  ♦ (t$*t/Vj-  rif)ar/V)V3t2«/VJ*  (W) 

M2/<r-U''$x  * ?./v  ?'$„  *?„;  * 2fx  ♦ >i%  *„  * 

The  boundary  condition  on  the  airfoil  surface  can  be  expressed  as 
v . Vs  * 0 
In  this  case 

t - i (i  ♦ tx)  r.  < fx;«  j fc,  *4-2  . 


Therefore,  on  the  aifrfoil  s uppepr  surface 
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r„ 


*Z  ! 2-0  • - s - !*»  ^(x  * vt)/v  j*€ 
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and 


'I-°  ’^'x-0  <fx-ef«xx.x-0 
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(B--5) 


The  discontinuity  in  y*  across  the  leading  edge  Mach  sheet  need  not 
be  evaluated  since  the  differential  equation  and  the  boundary  conditions 
on  the  airfoil  surface  are  sufficient  to  determine  its  value. 

The  first  order  boundary  condition  on  the  airfoil  s upper  surface 
for  a constant  vertical  acceleration  and  a steady  pitching  is  (from 
eq.  (8),  and  eq.  (8)  of  ref.  (5)) 


**|  * -®  - [ot  + (qx/^*  €1* 


2-0 

Comparing  eq.  (B-U)  and  (B-6)  it  can  be  seen  that  if  we  let 

4 

a • q - a ixo. 

Iu 


(B-6) 


(B-7) 


to  the  first  order  the  oscillating  motion  can  be  considered  as  the  sura 
of  a constant  vertical  accele -ation  a*Vi  a steady  pitching  motion. 
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The  second  order  boundary  condition  on  the  airfoil's  upper  surface 
for  a constant  vertical  acceleration  and  a steady  pitching  is  (from  eq. (9), 
and  eq.  (9)  of  ref.  (£)) 


where. 0^  denotes  the  first  order  solution  associated  with  pitching, 
and  0°  denotes  the  first  order  solution  associated  with  a constant 
vertical  acceleration.  Comparing  eqs.  (B-5)  and  (B  8)  it  can  be  seen 
that  the  second  order  boundary  condition  on  the  airfoil  surface  for  a 
slowly  oscillating  airfoil  can  be  considered  as  the  sum  of  a constant 
vertical  acceleration  arid  a steady  pitching  motion. 


The  second  order  partial  differential  equation  for  the  combined 
motion  of  accelerating  and  pitching  can  be  written  as 


2M  (x  ia 


vt  * **.t)  * *VV)(fcxx*$J 
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M-c  * f*  x *«/V  ; Y#z  /T' 


where  q and  a nave  been  replaced  by  a ia>  Eqs  (3-3)  and  (B-9) 
are  not  the  same;  they  agree  except  for  the  terms  2M  ? 

a 2M2  ./V  , a ixo  2M2  ztf^/V2,  and  a iu>  2M2  x <5  ./V2. 

d Zv  m XL  in  z v 

These  terms  do  not  contribute  to  the  part  of  the  second  order  solution 

which  is  cf  the  form  a iu  € h(x,y,z,t).  This  can  be  established  by 

analyzing  the  form  of  ftie  second  order  solution  in  the  same  manner  that 

the  form  of  the  accelerating  solution  was  analyzed  previously  in  this 

paper.  An  analysis  of  the  form  of  the  second  order  solution  for  the 

slowly  oscillating  airfoil  will  also  show  that  only  that  part  of  the 

second  order  solution  which  is  of  the  form  ico  €h(x,y,z,t)  will 

contribute  to  the  lifting  pressure  which  is  out  of  phase  with  the  in- 
stantaneous angle  of  attack,  and  Jicnce  to  the  aerodynamic  damping  of 
the  motion.  It  follows  that  to  the  second  order  the  damping  of  slowly 
oscillating  airfoils  is  directly  proportional  to  Cm  ♦ Cm  . 

& q 
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